X-ray views of neutron star low-mass X-ray binaries 
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A neutron star low-mass X-ray binary is a binary stellar system with a neutron star and a low- 
mass companion star rotating around each other. In this system the neutron star accretes mass 
from the companion, and as this matter falls into the deep potential well of the neutron star, the 
gravitational potential energy is released primarily in the X-ray wavelengths. Such a source was first 
discovered in X-rays in 1962, and this discovery formally gave birth to the "X-ray astronomy" . In the 
subsequent decades, our knowledge of these sources has increased enormously by the observations 
with several X-ray space missions. Here we give a brief overview of our current understanding of 
the X-ray observational aspects of these systems. 

PACS numbers: 



I. INTRODUCTION 

Cosmic objects emit in diverse electromagnetic wave- 
lengths from radio to 7-rays. Some of them (e.g., stars) 
are luminous primarily in a narrow wavelength range, 
while others (e.g., active galactic nuclei) emit in a broad 
energy spectrum. It is, therefore, essential to observe a 
celestial source in various wavelengths in order to fully 
understand its nature. Up to the middle of the twentieth 
century, visible light was the primary window for obser- 
vation. But in the last few decades, wavelength-based 
branches of astronomy have not only been developed, but 
also got matured. Moreover, many new types of sources 
have been discovered by the observations in wavelengths 
other than visible lights. One such discovery was made in 
X-rays using rocket-borne instruments in 1962 These 
instruments detected X-rays from sources outside the so- 
lar system for the first time, and specifically discovered a 
new type of source Scorpius X-l. This discovery brought 
a new era in astronomy, and formally gave birth to a new 
branch called "X-ray astronomy" . The new source Scor- 
pius X-l was the first observed X-ray binary, and more 
specifically it was the first detected low-mass X-ray bi- 
nary system containing a neutron star and a low-mass 
stellar companion. 

Unlike the optical (i.e., visible light) and radio detec- 
tion instruments, the X-ray instruments must be sent 
above the atmosphere for observations. This is because 
X-rays from space cannot reach the surface of the Earth. 
In 1960's, balloons and sub-orbital rockets were exten- 
sively used to lift the X-ray instruments. However, as 
these carriers are short-lived, satellites were thought to be 
ideal for X-ray observations. The first astronomy satel- 
lite Uhuru was launched by NASA in 1970. It identified 
over 300 discrete X-ray sources. NASA's Einstein satel- 
lite (launched in 1978) was the first fully imaging X-ray 
telescope put into space. Its few arcsecond angular reso- 
lution increased the sensitivity enormously, and made the 
imaging of extended objects, and the detection of faint 
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sources possible. The R OS A T satellite launched in 1990 
by the Germany/US/UK collaboration did the first X-ray 
all-sky survey using a highly sensitive imaging telescope, 
and detected more than 150000 objects. The scientific 
outcome of these key X-ray space missions, as well as 
the other past X-ray satellites, have revolutionized our 
knowledge of the X-ray sky. 

The current primary X-ray space missions are (1) 
NASA's Rossi X-ray Timing Explorer (RXTE); (2) Euro- 
pean Space Agency's XMM-Newton; (3) NASA's Chan- 
dra; and (4) Japan's Suzaku. They primarily operate 
in the w 1 - 10 keV range (i.e., soft X-rays), and con- 
tain a variety of instruments. The primary instrument of 
RXTE is a set of five proportional counters (called "Pro- 
portional Counter Array" or PCA) with large collecting 
area and very good time resolution (« 1 microsecond). 
Therefore, PCA is an ideal instrument to study the fast 
timing phenomena, as well as the continuum energy spec- 
trum. However, it is not an imaging instrument, and due 
to poor energy resolution, PCA is not very suitable to 
study spectral lines. This instrument operates in 2 — 60 
KeV band, although the effective collecting area becomes 
small above « 20 keV. Another very useful instrument of 
RXTE is "All-Sky Monitor" or ASM. It consists of three 
relatively low-sensitivity proportional counters with large 
field of view, and hence monitors the entire sky to detect 
transient X-ray sources. XMM-Newton contains X-ray 
telescopes, charge-coupled devices and high resolution 
"Reflection Grating Spectrometers (RGS)". Therefore, 
this satellite, that operates in » 0.2 — 12 keV range, is 
a good imager, and is ideal for spectral analysis. The 
Chandra observatory contains an X-ray telescope, charge- 
coupled devices, high resolution camera and transmission 
gratings. Its angular resolution (~ subarcsecond) is the 
best among all the past and current X-ray instruments, 
which makes it a unique imager. Chandra is also ideal 
to detect and study narrow spectral lines, and it oper- 
ates in w 0.1 — 10 keV. Finally, the Suzaku satellite also 
has X-ray telescopes and charge-coupled devices. Since 
it contains an additional "Hard X-ray Detector" (HXD), 
Suzaku can study the broadband energy spectrum. 

All these X-ray space missions, especially the current 
ones, have gathered a wealth of information about the 
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low-mass X-ray binary systems, and hence our under- 
standing of these sources has grown enormously. In this 
review, we will give a brief overview of the X-ray obser- 
vational aspects of neutron star low-mass X-ray binaries. 
We note that the reference list will not be exhaustive in 
this short review. 



II. NEUTRON STARS 

Neutron stars are extremely compact objects with a 
mass of about one solar mass (M ), but a radius of ~ 10 
km. The core density of neutron stars are 5 — 10 times the 
nuclear density, and their magnetic field can be anywhere 
between 10 — 10 15 Gauss. They can spin with a wide 
range of frequency values: very slow to several hundred 
Hertz. Study of these stars provides a unique opportu- 
nity to probe some extreme aspects of the universe, that 
cannot be achieved by other means. 

Although the neutron stars were discovered in 1967, 
their plausible existence was predicted in 1930's. Around 
this time, Subrahmanyan Chandrasekhar calculated the 
upper limit of the mass of a collapsed star (white dwarf) 
that is supported by the electron degeneracy pressure 
against the gravity. This limit is the celebrated "Chan- 
drasekhar mass limit", with a value of about 1.4M . If 
the mass of the collapsed star is greater than this limit, 
then the star is expected to collapse more. But if the 
mass is less than a certain value, the neutron degener- 
acy pressure will eventually balance the gravity to form 
a stable star, known as neutron star [2j,|3|. This remained 
as a theoretical calculation almost for three decades, af- 
ter which Hewish et al. |4J reported the discovery of a 
1.377 sec period radio pulsar at 81.5 MHz, and Gold [HQ 
showed that the pulsars are spinning neutron stars with 
large surface magnetic fields. 

The density of a neutron star increases by w 14 orders 
of magnitude from the surface to the centre. As a result, 
the constituents and their properties change dramatically 
with the radial distance (see, for example, Shapiro & 
Teukolsky Q). The dense core of a neutron star may con- 
tain pion condensates, deconfined quarks, or other exotic 
(and possibly unknown) forms of matter, which makes 
the study of this part of the star extremely interesting. 
Moreover according to a model, some of the neutron stars 
are actually "strange stars" made entirely of u, d and s 
quarks 8]. Unfortunately, in spite of many theoretical 
studies, the properties of matter in the core of a neutron 
star is not yet known, because such a dense degenerate 
matter cannot be created in the laboratory. Perhaps the 
only way to understand the nature of this high density 
matter is to constrain the theoretically proposed equation 
of state models of the stellar core. This can be done by 
measuring three independent structural parameters (e.g., 
mass, radius, spin frequency) of the same neutron star, 
because the stable stellar structure gives a single mass 
vs. radius curve for a known spin frequency and a given 
equation of state model [9[ . Measurement of the stellar 




FIG. 1: Sections in the orbital plane of the equipotential 
surfaces of a binary system with the mass ratio 2:1. The 
two gravitating bodies are rotating around each other, and 
the above mentioned potential includes the effects of both 
gravitational force and centrifugal force. The two parts of 
the figure-of-eight equipotential (containing the L\ point) are 
called Roche lobes, and the inner Lagrange point L\, having a 
saddle-like potential, allows the matter to pass from one lobe 
to the other (figure courtesy: Dipankar Bhattacharya). 

parameters is, therefore, a primary goal of the neutron 
star astronomy. Such a measurement, especially of the 
radius, is an extremely formidable task. This can be ap- 
preciated from the fact that we need to measure the ~ 10 
km radius of a neutron star with < 5% error, while the 
typical distance of the star is ~ 10 17 km. 

Neutron stars can emit in various wavelengths depend- 
ing on their properties and the systems they are in. They 
are observed as radio, X-ray and 7-ray pulsars (that may 
create pulsar wind nebulae), millisecond pulsars in bi- 
nary systems, magnetars, and so on, and can be found in 
supernova remnants, X-ray binaries, etc. In this review 
we will concentrate on neutron stars in X-ray binaries. 



III. LOW-MASS X-RAY BINARY 

A low-mass X-ray binary (LMXB) is an old binary 
stellar system (typical age ~ 10 9 years) with a low-mass 
companion star (< 1M Q ; main sequence star, evolved 
star or white dwarf) and a primary compact star (neutron 
star or black hole) rotating around each other. In such 
a system, the companion star fills its Roche lobe, i.e., 
the critical equipotential surface (see Fig. 1), and mat- 
ter from it flows towards the compact star. Because of 
the initial angular momentum, this matter cannot move 
radially towards the primary star, and follows slow spi- 
ral paths that form an accretion disk. In this review, we 
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A Low Mass X-Ray Binary: 4U 1820-30 




FIG. 2: Artist's impression of a neutron star LMXB 4U 1820- 
30 22]. This is a particularly small LMXB with the orbital 
period of 11 minutes (Sun and Earth are shown for an easy 
comparison of size). The companion is naturally a white 
dwarf star, because only such a star can have the size sim- 
ilar to that of companion's Roche lobe in these small systems 
(figure courtesy: Dany P. Page). 



consider the systems with neutron stars as the primary 
stars. In such neutron star LMXBs, the accreted matter 
eventually hits the stellar surface, and generates electro- 
magnetic radiation. The accretion disk also emits such 
radiation by viscous dissipation, and both these emis- 
sions are powered by the gravitational potential energy 
release. The inner part of the accretion disk, as well as 
the neutron star surface, emit X-rays, as their typical 
temperature is ~ 10 7 K. Fig. 2 gives an artist's impres- 
sion of such a system. In the next section, we will give 
the basic theory of a simple geometrically thin accretion 
disk in Newtonian gravity. Note that a more realistic 
approach may include complicated fluid dyanamical pro- 
cesses (e.g., Mukhopadhyay 2006 [10(; Mukhopadhyay et 
al. [ll|) and general relativistic effects (e.g., Mukhopad- 
hyay & Misra 2003 [HI])- Moreover, the magnetosphere of 
the neutron star can affect the accretion process depend- 
ing on the accretion rate and the magnetic field strength. 
Ghosh et al. [l3| and Ghosh & Lamb [HJ provided im- 
portant theoretical understanding of this process. Due to 
accretion induced angular momentum transfer, the neu- 
tron stars in LMXBs are expected to be spun up. This 
is consistent with the observed stellar spin frequencies 
(typically ~ 300 — 600 Hz). The magnetic fields of these 
neutron stars are relatively low (« 10 7 — 10 9 Gauss) com- 
pared to the expected fields of these stars when they are 
born, and several groups have worked on this problem 
of stellar magnetic field evolution [15l - l20l ] . Note that the 
low magnetic field allows the accretion disk to approach 
very close to the neutron star. This, in combination with 
the fact that X-rays can be detected both from the star 
and the inner part of the disk for many LMXBs, provides 
an excellent opportunity to study the extreme physics of 



and around neutron stars. However, the modeling of any 
observation of the neutron star vicinity should consider 
the special and general relativistic effects including light 
bending (e.g., Datta & Kapoor |21|). 

LMXBs are observed mostly in the disk, bulge and 
globular clusters of our galaxy. Some of them have also 
been discovered from the nearby galaxies, and the lo- 
cations much above our galactic disk plane. According 
to the catalogue of Liu et al. [Hj], 187 LMXBs (con- 
taining black hole or neutron star) were known from our 
galaxy, "Large Magellanic Cloud" and "Small Magellanic 
Cloud" up to that time. LMXBs emit primarily in X-ray 
wavelengths, and the optical to X-ray luminosity ratio 
is normally less than 0.1 for them. The central X-ray 
emission from these sources cannot be spatially resolved 
using the current X-ray imaging instruments, because the 
distances of even our galactic LMXBs are very large (typ- 
ically ~ 10 17 km). Therefore, in order to understand the 
neutron star LMXBs, we have to rely on their spectral 
and timing properties. In this review, we will describe 
some of these properties, as well as the information they 
provide. 



A. Theory of Accretion Disk 

In a neutron star LMXB, the companion mass is trans- 
ferred to the neutron star typically via a geometrically 
thin accretion disk. Therefore, in this section we will 
briefly discuss the basic theory of such a disk. Here we 
note that the standard model of steady thin accretion 
disks in Newtonian gravity was worked out by Shakura 
& Sunyaev |23| . 

In a thin accretion disk using cylindrical coordinates 
(r, 6, z), we expect that the 9 component vg should be 
the main component of velocity, and the z component 
v z = 0. The accreted matter should move around the 
neutron star with a nearly Keplerian speed (i.e., the an- 
gular speed Q » i/GM/r 3 ; M = mass of the neutron 
star), and a small radial inflow (speed = v r << Vg) 
caused by the effect of viscosity should exist. We con- 
sider d/d6 = because of the 0-symmetry of the system. 
Then the equation of continuity is 



§ + -|:(^) = o, 

ot r or 



(1) 



where E is the surface density of the accretion disk. Using 
this and the f9-component of the Navier-Stokes equation, 
we get 

jWo )+ i!_( Er 3 ni;r)=Ai (2) 

at r or 

where f2 = vg/r and A is the term involving viscosity. 
We note that 2Trr.dr.T,r 2 fl is the angular momentum of 
an annulus between r and dr. Hence the equation (2) 
multiplied with 27rr.dr gives how the angular momentum 
of the annulus changes by the viscous torque. Therefore, 

A = ^, (3) 
27rr dr J 
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where G(r) is the viscous torque at the radial distance 
r. It is normally assumed that the term dvg /dr gives the 
velocity shear. However for the accretion disk, 



dr dr dr 



(4) 



in which ft is associated with pure rotation, and rdil/dr 
gives the shear. Therefore, the viscous stress is /ird£l/dr, 
where p is the coefficient of viscosity. Then the viscous 
torque is 



G(r) = J rd9 J d 



zfj,r~—— = 2m/Er -r— 
dr dr 



(5) 



Here, v = fj,/p is the kinematic viscosity and p is the 
volume density. Therefore using the equations (2), (3) 
and (5), we get 

JWfi) + i|-(Sr 3 ^ r ) = (6) 
at r or r or dr 

The equations (1) and (6) govern the accretion disk dy- 
namics. If the disk is nearly Keplerian, then these two 
equations can be combined to get 



— = - — [r 1/2 — (i/Er 1/2 )] 
dt r dr dr 



(7) 



This equation gives the time evolution of the accretion 
disk. 

Let us now consider a special case, viz., a steady ( = 
0) geometrically thin accretion disk. From the equations 
(1) and (6), we get 



rT,v r = C\, 



Sr J rJw r - vY>r A — = C 2 , 
dr 



(8) 



(9) 



where C\ and C2 are constants. For a steady disk, the 
mass accretion rate is rh = —2irrHv r . Therefore from 
equation (8), 



m 

'2tt' 



(10) 



In order to calculate C2, we consider a boundary condi- 
tion 



™| =0 



(ii) 



Let us justify this condition for the following two cases: 
(1) when the disk is extended up to the stellar surface, 
and (2) when the star is well inside the disk, and the mat- 
ter falls from the disk inner edge towards the star quasi- 
radially. For a near-Keplerian disk, the angular speed 
increases with the decrease of r. For the first case, since 
the stellar angular speed is less than Keplerian (otherwise 
the star would break up) , the disk f2 must decrease with 
the decrease of r close to the star. Therefore, ^ must be 



zero at some point (assumed to be r = n n ) near the stel- 
lar surface. On the other hand, for the second case, the 
torque at the disk inner edge is expected to vanish (note: 
it may not be exactly true in reality) making 4P = 0. 
The region (between the star and the disk) in which the 
flow deviates considerably from Keplerian is called the 
boundary layer. From the equations (8), (9) and (11), 
we get 



2tt 



(GMr h 



a/2 



(12) 



where M is the stellar mass. Substituting this in the 
equation (9), and using the equations (8) and (10), we 
get, 



^ 771 r 

uT, = — 1 
3tt L 



(r-m/r) 



1/2, 



(13) 



This shows that the accretion rate is greater for larger 
viscosity. The viscous dissipation in the accretion disk 
causes the emission from the disk, and the source of this 
emitted energy is the gravitational potential energy re- 
lease. This viscous dissipation rate per unit volume of the 
accretion disk is pr 2 (dfl / dr) 2 . Integrating this over the 
disk thickness and assuming Keplerian angular speed, we 
get the rate of energy emitted per unit area of the disk 
from equation (13): 

dE 3GMm 1/2 

-^^Hi-^An. (w) 

Therefore, the total power emitted from a disk extending 
to infinity is 

GMm 



from r 



L = 



2ru 



(15) 



This shows that, for a steady geometrically thin accre- 
tion disk in Newtonian gravity, half of the gravitational 
potential energy is emitted by viscous dissipation in the 
disk, and the rest half remains stored in the accreted 
matter as the kinetic energy, and gets dissipated in the 
boundary layer. 

It is expected that a geometrically thin Keplerian disk 
should emit like a blackbody. Therefore, the temperature 
T Q g(r) of the disk surface can be calculated from equation 
(14): 

Teff (r)={^p[l-(r in /r)^]}V4 ! (16) 
87rr J er 

where a is the Stefan-Boltzmann constant. As the black- 
body temperature varies with the radial distance, the 
spectrum from a disk is known as multicolor blackbody. 
For luminous LMXBs, the X-rays from the neutron star 
boundary layer and the inner part of the accretion disk 
can substantially irradiate the disk surface. This irradi- 
ation temperature Ti rr (r) (which is blackbody tempera- 
ture) is proportional to r -1 / 2 [24]]. Therefore, the net ef- 
fective disk temperature is T disk (r) = [r e 4 ff (r)+7]4 r ( r )]i/4 ; 
in which T e g(r) dominates for the inner portion (that 
emits X-rays) of the disk, and T irr (r) dominates far away 
from the neutron star. Note that the latter can have 
significant consequences for the disk instability [24|. 
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FIG. 3: Long-term X-ray lightcurve (i.e., observed intensity 
profile) of a transient neut ron star LMXB 1A 1744-361 show- 
ing its two outbursts [l!2l |. The intensity is expressed in the 
photon count rate detected with the PCA instrument of the 
RXTE satellite. 



IV. SELECTED FEATURES 

A. Outbursts 

A class of neutron star LMXBs are called transients. 
These sources have two distinctly different intensity 
states: (1) the quiescent state in which the source re- 
mains for months to years; and (2) the outburst state 
which continues for weeks to months (sometimes for 
years; e.g., for the source KS 1731-260). The luminos- 
ity of the source typically increases by several orders of 
magnitude as it evolves from the quiescent state into an 
outburst. Fig. 3 shows the long term intensity profile of 
such a transient 1A 1744-361, exhibiting two clear out- 
bursts. The transient nature of these LMXBs is normally 
attributed to an accretion disk instability, that causes 
high accretion rates for certain periods of time, and al- 
most no accretion during other times (see King 2001 25 1 
for a review). These sources are often not detectable in 
their quiescent states, and in some cases when they are 
detectable, the poor statistics due to the small number 
of observed photons hinders a detailed study. For most 
of the scientific purposes, it is therefore preferable to ob- 
serve these transients during their outbursts. New tran- 
sients are also typically discovered during the outburst 
states. However, it cannot normally be predicted when a 
transient will come out of its quiescent state. Therefore, 
instruments with large field of view that can monitor the 
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FIG. 4: X-ray intensity profile of a typical thermonuclear X- 
ray burst observed from an LMXB. After the onset of the 
burst, the X-ray intensity increases quickly by more than an 
order of magnitude (compared to the non-burst intensity cor- 
responding to the gravitational energy release), and then de- 
cays relatively slowly as the neutron star surface cools down. 



entire sky almost continuously are required to find the 
transients in their outbursts. The RXTE ASM is such 
a very useful instrument, that has discovered many new 
transient LMXBs, as well as detected subsequent out- 
bursts of them. 

However, there are scientific reasons that motivate the 
observations of transient neutron star LMXBs in their 
quiescent states. For example, during the quiescence, 
the neutron star surface should be the primary X-ray 
emitting component, which can be modelled to constrain 
the stellar radius. Besides, a long outburst heats up the 
neutron star crust taking it out of thermal equilibrium 
with the rest of the star. Therefore, tracking the LMXB 
system going into the quiescence, and the analysis of the 
quiescent X-ray data can provide a unique opportunity to 
probe some aspects of the neutron star extreme physics 
(e.g., Cackett et al. [Hj]). For example, such a study 
allows one to compare the observed neutron star cool- 
ing rate with the theoretical cooling curve, that involves 
the physics of stellar constituents. Note that the satel- 
lites (e.g., Chandra, XMM-Newton) with high sensitivity 
imaging instruments are used to observe the sources in 
quiescence. 



B. Thermonuclear X-ray Bursts 

Many neutron star LMXB systems show a spectacu- 
lar eruption in X-rays every few hours to days. Dur- 
ing these bursts, the observed X-ray intensity goes up 
sharply in ss 0.5 — 5 seconds, and decays relatively slowly 
in ss 10 — 100 seconds (see Fig. 4; Strohmayer & Bild- 
sten The typical energy emitted in a few seconds 

during such a burst is ~ 10 39 ergs (note that the Sun 
takes more than a week to release this energy). These 
bursts are known as type-I or thermonuclear bursts, and 
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were first discovered in 1970's [28|, 129J. Soon after their 
discovery it was realized that they originate from recur- 
rent unstable nuclear burnin g of accreted matter on the 



the subsequent j3 decay lifetimes. Hence, the hydrogen 
burns via the "hot" CNO cycle of Fowler & Hoyle j35j: 



neutron star surfaces [2JJ, [30|, |3l| . The neutron star sur- 
face origin was supported by the observational fact that 
the burst emission area (estimated from the energy spec- 
trum of the bursts) matched with the expected surface 
area of a neutron star [32| . Although, a competetive 
emission area could be the inner portion of the accretion 
disk, nuclear ignition in order to produce the observed 
bursts cannot occur in such a rare medium. But what 
shows that the nuclear burning causes the bursts? Be- 
fore coming to that point, we would like to convince the 
readers that such a nuclear burning has to be unstable. 
This is because, while the nuclear energy released per 
nucleon by fusion is only a few MeV, the gravitational 
energy release by a nucleon falling on the neutron star 
surface from a large distance is GMm nuc \ con / R (M and 
R are neutron star mass and radius respectively), which 
is typically ~ 200 MeV. Therefore, the energy released 
from a stable burning of the accreted matter would be 
entirely lost in the observed gravitational energy. So the 
nuclear burning can power the bursts, only if the ac- 
creted matter accumulates on the neutron star surface 
for a long time, and then burns quickly in an unstable 
manner. This instability appears if the nuclear energy 
generation rate is more temperature sensitive than the 
radiative cooling rate. For an observed data set from 
a given source, the ratio of the total energy content in 
all the bursts to that in the non-burst emission roughly 
tallies with the expected nuclear energy release (per nu- 
cleon) to the expected gravitational energy release (per 
nucleon) ratio (33| This is one of the strongest evidences 
of the nuclear energy origin of the bursts. Here we note 
that, only a year before these bursts were discovered, 
Hansen and Van Horn [34[ predicted that thin hydrogen 
and helium layers on neutron star surfaces could be sus- 
ceptible to the above mentioned instability 27 1. 
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The accreted material consists of mostly hydrogen, 
some helium and a small amount of heavier elements, 
except when the donor companion is a white dwarf. This 
matter, that accumulates on the neutron star surface, un- 
dergoes hydrostatic compression as more and more ma- 
terial piles on, and the ignition density and temperature 
are reached in a few hours to days. The burning typi- 
cally happens at a depth of 10 m at a column density 
of ~ 10 8 gm cm" 2 . The compression rate and the igni- 
tion depends on the accretion rate, which sets four dis- 
tinct regimes of burning [33j . Regime 1: in most cases, 
the temperature of the accumulated material exceeds 10 7 
K, and hence the hydrogen burns via the CNO cycle in- 
stead of the proton-proton (pp) cycle. This causes mixed 
hydrogen and helium bursts triggered by the unstable hy- 
drogen ignition, which happens for the accretion rate per 
unit neutron star surface area rh < 900 gm cm~ 2 sec -1 
(^cno/0-01) 1 / 2 [1^. Here Zcno is the mass fraction of 
CNO. Regime 2: at higher temperatures (T > 8 x 10 7 
K), the proton capture timescale becomes shorter than 



C(p, 7 ) 13 N(p, 7 ) 14 0(/3+) 14 N(p, 7 ) 15 0(/3+) 15 N(p, a) 12 C. 



In this case, the hydrogen burning happens in a ther- 
mally stable manner (i.e., without triggering a ther- 
monuclear burst) simultaneously with the accumulation 
of matter. This allows a helium layer to build up be- 
low the hydrogen layer. For the rh range 900 gm cm 
sec" 1 (ZcNo/0.01) 1 / 2 < m < 2 x 10 3 gm cm" 2 sec" 1 
(ZcNo/0-01) 13 / 18 , the hydrogen becomes entirely burned 
before this helium is ignited [13, HH . So when the he- 
lium ignition happens, a short (~ 10 sec) but very in- 
tense burst occurs by the unstable triple-alpha reaction 
of the pure helium (3a — > 12 C). Such a burst is called 
a "helium burst" . Many of these bursts are so intense, 
that the local X-ray luminosity in the neutron star at- 
mosphere may exceed the Eddington limit (for which the 
radiative pressure force balances the gravitational force), 
and the photospheric layers may be lifted off the neutron 
star surface. Such bursts are called photospheric radius 
expansion (PRE) bursts. Regime 3: for rh > 2 x 10 3 gm 
cnr 2 sec -1 (Zcno/0.01) 13 / 18 , hydrogen burns via the 
"hot" CNO cycle in a stable manner (as for the regime 
2), but enough unburnt hydrogen remains present at the 
time of helium ignition. This is because at this higher ac- 
cretion rate, the helium ignition conditions are satisfied 
much sooner. Therefore in this regime, mixed hydrogen 
and helium bursts are triggered by the helium ignition. 
During such a burst, the thermal instability can produce 
elements beyond the iron group [36l - l40| via the rapid- 
proton (rp) process of Wallace & Woosley jUl . This pro- 
cess burns hydrogen via successive proton capture and /3 
decays. The long series of (3 decays makes these bursts 
typically much longer (~ 100 sec) than the helium bursts. 
Regime 4'- finally at a very high accretion rate (compara- 
ble to the Eddington limit), the helium burning temper- 
ature sensitivity becomes weaker than the cooling rate's 
sensitivity [IH, HI]- Therefore in this regime, the sta- 
ble burning sets in, and the thermonuclear bursts do not 
occur. 

Apart from type-I X-ray bursts, another type of ther- 
monuclear X-ray bursts with much larger recurrence time 
(~ years) has been observed. These are called super- 
bursts, and they differ from the type-I bursts also in the 
decay time (~ 1 — 3 hours) and the amount of released 
energy (~ 10 42 ergs). The first superburst was discovered 
by Cornelisse et al. fiij from the known type-I burster 4U 
1735-44. So far about 15 superbursts have been detected 
from 10 sources [45| . It is believed that these bursts of 
large ffuence happen by the 12 C + 12 C fusion reaction at 
a column depth of 10 12 gm cm" 2 [45T - l47l | . However, 
currently our understanding of superbursts is not as clear 
as that of type-I bursts. 

Thermonuclear bursts provide a unique opportunity 
to study some aspects of extreme physics. Moreover, 
such a study brings together several branches of physics: 
nuclear physics, astrophysics, general theory of relativ- 
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ity, magnetohydrodynamics, etc. These bursts can also 
be useful to constrain the neutron star parameters, and 
hence to understand the nature of the stellar core mat- 
ter at supranuclear densities (see the spectral and timing 
properties of these bursts in later sections). Therefore, 
the thermonuclear bursts belong to a truly multidisci- 
plinary field, which is not yet well understood. One such 
poorly explored aspect of these bursts is the thermonu- 
clear flame spreading on the neutron star surfaces. Burst 
is ignited at a certain location on the stellar surface, and 
then the burning region expands to cover the entire sur- 
face [3 |4i| . The theoretical modeling of thermonuclear 
flame spreading including all the main physical effects is 
extremely difficult, and only recently Spitkovsky et al. 
f40l | have provided some insights in this field. Although 
these authors have ignored several physical effects (mag- 
netic field, strong gravity, etc.) in their model, they have 
considered the effects of the Coriolis force, which is im- 
portant as the bursting neutron stars are rapidly spin- 
ning (spin frequency ~ 300 — 600 Hz). The observational 
study of flame spreading is also no less difficult than the 
theoretical study. This is because, this spreading mostly 
happens within « 1 sec during the burst rise, and in or- 
der to study it observationally, one needs to detect and 
measure the evolution of the burst spectral and timing 
properties during this short time. Only recently this has 
been possible for a few bursts by analyzing the RXTE 
PC A data [50l - l52j . Moreover, an indication of the effects 
of the Coriolis force on the flame spreading (as predicted 
by Spitkovsky et al. [Hj]) has also been found [53( | . 
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FIG. 5: X-ray intensity dips observed from a neutron star 
LMXB 1A 1744-361 )112]. The separation between two sets 
of consecutive dips gives the orbital period of the binary sys- 
tem (see text). Because of the gap in the observed data, one 
cannot be sure whether the two sets of dips shown in this 
figure are consecutive. Hence, these two sets of dips may be 
the orbital period or a multiple of that. 



D. Dips and Eclipses 



C. Type-II Bursts 

A series of bursts, known as type-II X-ray bursts, are 
observed from neutron star LMXBs in quick succession. 
To the best of our knowledge, only two such systems 
(rap id burster and GRO J1744-28) exhibit these bursts 
[5J-|56j. The intervals between the bursts vary between 
seconds to ~ 1 hour, and the burst duration is in the 
range of seconds to minutes. The typical emitted energy 
range during a burst from the rapid burster is ~ 1 x 10 38 
7 x 10 40 



ergs to 



ergs 



571. These bursts often show 



successive peaks and flat tops. Note that a source is not 
always active for type-II X-ray bursts, and a series of such 
bursts are observed in certain periods of time. type-II X- 
ray bursts are distinctly different from the type-I bursts 
in many ways. For example, The latter bursts show a 
spectral softening during burst decay, while the former 
bursts lack this signature. type-II bursts are believed to 
be caused by spasmodic accretion, possibly due to an ac- 
cretion instability [57|- Therefore, these bursts are pow- 
ered by the gravitational energy, while the type-I bursts 
are powered by the nuclear energy. The recent models 
of the type-II bursts include the effects of magnetosphere 
58] and intermittent advection-dominated accretion flow 
59| . However, the origin of the accretion instability is not 
yet well understood. 



Some neutron star LMXBs show periodic temporary 
decrease in X-ray intensity. This timing feature can be 
divided into two distinct categories: eclipses and dips. 
Both of them occur with the binary orbital period, and 
hence provide an excellent way to measure this period 
[60j |. Partial and total eclipses are believed to be caused 
by the obscuration of the X-ray source by the companion 
Such an obscuration requires an observer's line- 



61 



star 

of-sight close to the plane of the binary system (i.e., close 
to the accretion disk plane). Therefore, the observation 
of eclipses from a source implies that its inclination angle 
is large (> 75°; Frank et al. [62|). This inclination an- 
gle i is measured from the direction perpendicular to the 
accretion disk plane. The dips are periodic but irregu- 
lar, and typically are more pronounced at lower energies 
|6l| . The shape, depth and duration of each dipping in- 
terval widely vary from cycle to cycle, and within such 
an interval the observed intensity varies dramatically on 
a timescale of a few seconds (see Fig. 5). These dips 
are believed to be caused by a bulge on the edge of the 
accretion disk at the region where the gas stream from 
the companion star hits the disk (see Fig. 2 for a view 
of this impact region; White & Swank [63| ). Therefore, 
dips imply i > 60° 62j . The detailed nature of this bulge 
is still a matter of debate, and spectral analysis can shed 
light on this. 
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FIG. 6: Broad relativistic iron (Fe Kq) spectral line from the 
inner accretion disk of a neutron star LMXB Serpens X-l [8(|. 
The x-axis shows the X-ray energy, and the y-axis gives the 
observed intensity in excess to the best-fit continuum spec- 
tral model. The data points (shown with error bars) clearly 
shows an asymmetric spectral line, while the dotted profile is 
a relativistic model of spectral line that fits the observed line 
well. 



V. SPECTRAL ASPECTS 
A. Continuum Spectrum 

1. Persistent Emission 

The study of the continuum spectrum of the persistent 
emission (i.e., non-burst, non-dipping, non-eclipsing) 
from the neutron star low-mass X-ray binary systems can 
be very useful to understand the various components of 
the accretion flow and their structures, including the flow 
in the strong gravity region near the neutron star. Such 
a spectrum can be typically fitted with two spectral com- 
Mitsuda et al. 



ponents. Mitsuda et al. |64j used a single-temperature 
blackbody plus a multicolor blackbody model. The for- 
mer component is expected to be produced from a bound- 
ary layer between the accretion disk and the neutron 
star, while the latter may be produced from the accretion 
disk. In an alternative two-component model, White et 
al. [H, |6(| used a Comptonization component instead of 
a multicolor blackbody. They assumed that the Comp- 
tonization component was originally the disk component 
which was dominated by the Comptonization processes. 
In the later years, it has been found that a Comptoniza- 
tion component is often required to fit the observed per- 
sistent emission continuum spectrum, and an additional 
blackbody is needed when the source is more luminous 
and spectrally softer (see, for example, Bhattacharyya et 
al. |67j). Apart from these, there is evidence of an ex- 
tended "accretion disk corona (ADC)" for some sources, 



that can be modelled with a power law [68j. However, 
we note that much work needs to be done in order to 
build a self-consistent spectral model involving all the X- 
ray emitting and absorbing components, and to fit the 
observed spectrum with such a single model. 



2. Thermonuclear Burst Emission 

Spectrum of a thermonuclear X-ray burst is obtained 
usually by subtracting the preburst emission spectrum 
(assumed to be powered by the gravitational energy re- 
lease) from the observed spectrum during the burst [69j j . 
The continuum burst spectra are thermal, and in most 
cases they are quite well described with a blackbody 
model [12, [z3| • The temperature of this blackbody natu- 
rally increases rapidly during the burst intensity rise, and 
then decreases relatively slowly during the burst decay. 
If the entire neutron star surface emits like a blackbody 
of the same temperature at a given time, then the neu- 
tron star radius (i? bs) can, in principle, be measured 
from the following formula: i? bs = d[F \, s / {oT^)] 1 / 2 . 
Here, F b s and T b s are the observed flux and temper- 
ature respectively, and d is the source distance, which 
may be known by an independent technique. Note that 
the measurement of the neutron star radius can be use- 
ful to understand the nature of the dense matter of the 
stellar core. However, in order to do such measure- 
ment, one needs to include the following two effects. (1) 
The electron (Compton) scattering in the neutron star 
atmosphere shifts the blackbody spectrum (originated 
from the burning layer) towards the higher energies, and 
changes its shape slightly [7lM75| . Such a blackbody is 
called a "diluted blackbody", and its temperature (T co \, 
i.e., the color temperature) is related with the blackbody 
temperature (Tbb) by the color factor /. (2) The surface 
gravitational redshift shifts the temperature towards a 
lower value. Therefore, the observed temperature T Q b s is 
related to Tbb by Tbb = T Q i, s (l + z) / f . Correspondingly, 
the actual neutron star radius -Rbb is then related to the 
observed radius R b s by -Rbb = -Robs/ 2 /(l + z). There- 
fore, in order to measure the neutron star radius from 
the continuum burst spectra, one needs to independently 
estimate or measure the color factor / and the surface 
gravitational redshift factor 1 + z. The former can be 
estimated from the temperature, surface gravity and the 
atmospheric chemical composition [76| . while the latter 
can be calculated from the neutron star radius-to-mass 
ratio. This ratio can be independently measured from 
the surface atomic spectral lines or the burst oscillations 
(see later). 

The evolution of the continuum spectrum of a ther- 
monuclear burst can be best measured using the PCA 
instrument of the RXTE satellite. Such a measurement 
during the burst rise will be extremely useful to study the 
thermonuclear flame spreading. The chemical composi- 
tion of the neutron star atmosphere may change during 
the burst. Therefore, and since the burst temperature 
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also evolves, the observed radius i? bs changes substan- 
tially, and apparently irregularly during the burst decay 
69] . An understanding of these changes will be essential 
to (1) measure the neutron star radius using the burst 
continuum spectrum, (2) understand the nuclear physics 
of bursts, and (3) probe the magnetohydrodynamics of 
the burning layer and the atmosphere. 



3. Dipping Emission 

In order to probe the nature of the bulge that causes 
dips, it is essential to understand how the continuum 
spectrum changes as the source goes into the dipping 
interval. A simple increase in the neutral absorber col- 
umn density from its persistent emission value cannot 
explain the dip spectra. Currently, there are two com- 
peting spectral models to explain dips. (1) The "pro- 
gressive covering" model (e.g., Church et al. [68[): in 
this approach, the X-ray emission is considered to orig- 
inate from a point-like blackbody and a power-law from 
an extended component (e.g., a corona). The dipping is 
assumed to be caused by the partial and progressive cov- 
ering of the extended component by a neutral absorber. 
(2) The ionized absorber model (e.g., Boirin et al. [z3): 
in this approach, the dips are caused by the increase in 
column density and the decrease in ionization state of an 
ionized absorber, along with a neutral absorber column 
density increase. This model can naturally explain some 
of the observed narrow spectral features. 



B. Line Spectrum 

1. Lines from High Inclination Sources 

Many of the high inclination neutron star LMXBs (for 
example, dipping sources) exhibit narrow spectral ab- 
sorption features (e.g., Boirin et al. 78J; Parmar et al. 
[79j | : Diaz Trigo et al. 1801 ). as well as narrow recombi- 



nation emission lines [81| . The most prominent of these 
lines are Fe XXV and XXVI absorption lines. This spec- 
tral features are believed to originate from the material 
above the accretion disk, and can be observed only when 
the observer's line-of-sight passes through this material. 
This is possible only for the high inclination sources. 
The narrow spectral features, therefore, can be very use- 
ful to constrain the structures and properties (e.g., ion- 
ization states, relative abundances of various elements) 
of the photoionized plasma above the accretion disks in 
these systems [8l[. Note that, although the spectral fea- 
tures are mostly seen from high inclination LMXBs, this 
plasma must be common to all LMXBs, and hence these 
features will be useful in understanding LMXBs in gen- 
eral. 



2. Line from Disk 

Broad asymmetric spectral iron emission line has been 
observed from many accreting supermassive and stellar- 
mass black hole systems [HI, [83[ . It is believed to origi- 
nate from the inner portion of the accretion disk, when 
the disk is irradiated by a hard X-ray source (e.g., an 
accretion disk corona). Such an irradiation creates a "re- 
flection" component of the X-ray spectrum. When the 
hard X-ray photon enters the colder disk, it is either scat- 
tered out of the disk, or destroyed by Auger de-excitation, 
or reprocessed into a fluorescent line photon that escapes 
the disk |84J. The observed iron Ka line near 6.4 keV 
is the strongest fluorescent line from such a system. As 
the matter in the inner accretion disk rotates very fast, 
and is close to the black hole, Doppler effect broadens 
the line, special relativistic beaming makes it asymmet- 
ric, and the gravitational redshift shifts it towards the 
lower energies. Therefore, these lines can be useful to 
study the accretion flow close to the black hole, to test 
the theory of general relativity, and to measure the black 
hole spin. 

For long time it has been known that the persistent 
spectra of many neutron star LMXBs contain broad iron 
emission line [851 ] . It was suspected that such a line might 
have an origin similar to that of the broad iron line from 
black hole systems. However, the asymmetry of these 
lines was not established, and hence their origin was not 
well understood. Recently, Bhattacharyya & Strohmayer 
[86| have, for the first time, established the inner disk ori- 
gin of such a line from the neutron star LMXB Serpens 
X-l (see Fig. 6). This was later supported for several 
other sources (e.g;, Cackett et al. [87j ; Pandel et al. [88[ ; 
Di Salvo et al. [89[). Therefore, the high precision ob- 
servations of these broad iron Ka line from neutron star 
LMXBs could be useful to constrain the stellar parame- 
ters (radius, angular momentum) and to probe the inner 
accretion disk. However, the evidence of the "reflection" 
component of the X-ray spectrum is still marginal [89| . 
which the theory and the future observations should ad- 
dress. 



3. Line from Neutron Star 

Observing atomic spectral lines from the surface of a 
neutron star has been one of the long-standing goals of 
the neutron star astrophysics. This is because, such an 
observation and the identification of the line provide the 
cleanest way to measure the neutron star radius-to-mass 
ratio using the following formula (for the Schwartzschild 
spacetime appropriate for a non-spinning neutron star): 
E cm /E ohs = l+z= [1-{2GM/Rc 2 )}- 1 ' 2 [Il|9l|. Here, 
E cm is the emitted energy of the line photons, _E bs is 
the observed energy of these photons (different from E cm 
because of the surface gravitational redshift), and M and 
R are the mass and the radius of the neutron star respec- 
tively. Note that the measurement of the stellar param- 
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eter R/M can be useful to understand the nature of the 
dense core matter of neutron stars. As the neutron stars 
in LMXB systems typically spin with high frequencies 
(~ 300 — 600 Hz), Dopplcr effect and special relativistic 
beaming should make the lines broad and asymmetric 
(similar to the lines from disks mentioned in the previ- 
ous section). Measuring R/M accurately from such a 
broad line (for which the definition of £" bs is not cer- 
tain) using the above formula is difficult. Bhattacharyya 
et al. [9l| have recently suggested a way to measure the 
stellar R/M using the above formula from such lines. 
These broad asymmetric lines can have low-energy and 
high-energy pea ks ( or dips for absorption lines). Bhat- 
tacharyya et al. [91| have found that the ratio of energy 
contents of these peaks (or, dips) contains information 
about the general relativistic frame-dragging. Therefore, 
surface atomic lines can be used to test the theory of 
general relativity in the strong gravity regime. 

Surface lines are likely to be detected from thermonu- 
clear X-ray bursts and hence from neutron stars in LMXB 
systems. This is because: (1) the continuous accretion 
and the radiative pressure may keep the line-forming 
heavy elements in the stellar atmosphere for the time- 
duration required for the line detection; (2) the relatively 
low magnetic fields (~ 10 7 — 10 9 Gauss) of the neutron 
stars in LMXBs should make the line identification easier 
by keeping the magnetic splitting negligible; (3) the lower 
magnetic field does not complicate the modeling of the 
atmosphere; (4) the luminous bursts give good signal-to- 
noise ratio; and (5) during the bursts, typically 90% of 
the total emission originate from the neutron star sur- 
face, and hence the surface spectral line may be largely 
free from the uncertainty due to the other X-ray emission 
components, such as the accretion disk. 

Indeed, the only plausible observation of the neutron 
star surface atomic spectral lines was from an LMXB 
EXO 0748-676 [92|. These authors detected two signif- 
icant absorption features in the XMM-Newton RGS en- 
ergy spectra of the thermonuclear X-ray bursts, which 
they identified as Fe XXVI and Fe XXV n = 2 - 3 lines 
with a redshift of z = 0.35 (see Fig. 7). In order to de- 
tect these lines, they had to combine the spectra of 28 
bursts. The identification of these two lines as the surface 
atomic spectral lines is reasonable for the following rea- 
sons: (1) the early burst phases (when the temperature 
was higher) showed the Fe XXVI line, while the late burst 
phases (when the temperature was lower) showed the Fc 
XXV line. This is consistent with qualitative tempera- 
ture dependence of the iron ionization states. (2) The 
narrow lines were consistent with the relatively slow spin 
rate of the neutron star (45 Hz; Villarreal & Strohmayer 
[93|]). (3) Both the lines showed the same surface gravi- 
tational redshift, as expected. This redshift is reasonable 
for the surface of a typical neutron star. The absorp- 
tion of some of the bursts photons (originated in deep 
burning layers) by the iron ions in the upper atmosphere 
might give rise to these redshifted lines. Chang et al. [§4[ 
showed that the observed strength of the iron lines could 



be produced by a neutron star photospheric metallicity, 
which was 2 — 3 times larger than the solar metallicity. 
However, these lines were not significantly detected from 
the new data from the later observations of the same 
source. Although, this casts doubts in the reality of these 
lines, it may be caused by the change in the photospheric 
conditions that might weaken the lines. Therefore, the 
evidence of the existence of surface atomic lines currently 
remains uncertain, which, given the importance of these 
lines, provides the motivation for the future generation 
X-ray instruments. 



VI. TIMING ASPECTS 

A. Color-color diagram 

A convenient way to study how the persistent spectrum 
of a neutron star LMXB changes with time, and how 
it is correlated with various timing properties and the 
occurrence of other features (e.g., thermonuclear X-ray 
bursts), is to compute the color-color diagram. In order 
to do this, typically five photon energies with increasing 
values are considered: hv\, hv2, hfs, hv± and hv^. Let 
us assume that the observed background-subtracted pho- 
ton count numbers in the ranges hv\ — h^, hv-i — hvj,, 
hvz — hv± and hv^ — hv§ are iVi , N2 , N3 and N± respec- 
tively. Then we define soft color = N2/N1 and hard color 
= N4 /N3 . The color-color diagram is a plot of hard color 
vs. soft color, which shows the change of spectral hard- 
ness of the source. Neutron star LMXBs can be divided 
into two classes based on their color-color diagrams: the 
so called "Z sources", that show distinct Z-like curves 
(Fig. 8), and the "atoll sources", that exhibit C-like 
curves (Fig. 9). These two types of curves are not only 
different by shape, but also they are correlated with dis- 
tinctly different timing properties [95[ . Apart from these 
differences, the Z sources are very luminous (luminosi- 
ties are close to Eddington), while the luminosities of the 
atoll sources typically remain between 0.01 — 0.2 of the 
Eddington luminosity (95|. The three branches of a Z 
curve, from upper to lower, are called horizontal branch 
(HB), normal branch (NB) and flaring branch (FB). The 
lower part of a typical atoll curve is called the banana 
state, while the upper part is known as the island state. 
The X-ray luminosity typically increases from HB to FB 
for a Z source, and from island state to banana state for 
an atoll source. We note that, although the various char- 
acteristics of the color-color diagrams are observationally 
very robust, most of them are not yet theoretically well 
understood. Accretion rate is known to be a parame- 
ter that changes along the Z curves and the atoll curves. 
But, given the complex properties of these curves, more 
parameters must be involved in the fundamental level. 
It is also not understood what, apart from the accretion 
rate, distinguishes between the Z and atoll properties. 
The recent discovery of the first transient Z source, that 
shows substantially changed properties at lower luminosi- 
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ties, holds the promise to take our understanding of these 
sources to a higher level [§6j]. 

B. Regular Pulsations 

Most of the radio pulsars (that are spinning neutron 
stars emitting periodic radio pulses) have spin period 
> 10 milliseconds, and high magnetic field (typically 
~ 10 12 — 10 13 Gauss). These are relatively young neu- 
tron stars. This is because, as a pulsar spins down (due 
to magnetic dipole radiation) and its magnetic field de- 
cays (in timescales of several million years) , it eventually 
ceases to be a pulsar [93]- However, a small fraction of 
the radio pulsar population is rapidly spinning (period 
< 10 ms), old (age ~ billions of years), and having weak 
surface magnetic fields (10 s — 10 9 Gauss). Since many 
of these pulsars are in binary systems, it has been be- 
lieved for long time that these neutron stars were spun 
up by the accretion-induced angular momentum transfer 
in LMXB systems [97j |. This model would be supported 
by the observation of millisecond-period X-ray pulsations 
from neutron star LMXBs, as such pulsations would be 
a missing link between the millisecond radio pulsars and 
LMXBs. But years of searching did not yield the detec- 
tion of any such coherent millisecond-period X-ray pul- 
sations. Finally in 1998, the first millisecond X-ray pul- 
sations (with the frequency « 401 Hz) were discovered 
from the transient LMXB SAX J1808.4-3658 with the 
RXTE satellite [98|. As X-ray pulsations (like the ra- 
dio pulsations) measure the neutron star spin frequency, 
this showed that the neutron star in this source spins 
with s=s 401 Hz, and hence, indeed the neutron stars in 
the LMXBs can be spun up. Several such sources have 
been discovered after SAX J1808. 4-3658 (see, for exam- 
ple, Wijnands 2006 [99|; Watts et al. 2009 (lQOj . Fig. 
10 shows the power spectrum of one such LMXB. The 
narrow peak at the frequency of 442.36 Hz shows the 
pulsations. Note that the computation of a power spec- 
trum is a convenient way to detect timing features. Such 
a spectrum (power vs. frequency) is a Fourier transform 
of the intensity profile (photon counts vs. time). For ex- 
ample, if Xj(j = 0,...,N — 1) is the number of photons 
detected in bin j, then the discrete Fourier transform 
Pk{k — —N/2, N/2 — 1) decomposes the intensity into 
N sine waves. The direct and the inverse Fourier trans- 
forms arc then given by: 

N-l N/2-1 
3=0 k=-N/2 

In the LMXBs containing these pulsars, the accretion 
flow is channeled by the neutron star magnetic field to 
the magnetic polar cap. The X-ray hot spot on the neu- 
tron star created by this channeled flow rotates around 
the stellar spin axis (if the magnetic axis and the spin 
axis are misaligned) and gives rise to the observed X- 
ray pulses. Since the gravitational energy release from 



the accreted matter powers the X-ray pulsations, these 
sources are called "accretion-powered millisecond pulsars 
(AMPs)" . Since the first discovery of these pulsars, it in- 
trigues the astrophysicists how these sources are different 
from other LMXBs, and why all the LMXBs do not ex- 
hibit such regular pulsations. The recent discovery of 
several intermittent AMPs (that show sporadic regular 
pulsa tions) holds the promise to resolve these problems 
[l00j |. Better understanding of the AMPs will be very 
important to probe the evolution of LMXBs, as well as 
to model the magnetically channeled flows in the strong 
gravity regions near the neutron stars. Modelling of this 
flow, its energy spectrum and the timing properties will 
be useful to constrain the neutron star par amet ers, and 
hence the nature of the stellar core matter llOll . 



C. Burst Oscillations 

High frequency timing features are observed from 
many (so far about 20) neutron star LMXBs during ther- 
monuclear X-ray bursts [13] . They are called "burst 
oscillations", and were first discovered with the RXTE 
PC A instrument in bursts from the LMXB 4U 1728-34 
[l02j | . It was soon realized that this timing feature might 
be originated from an aymmetric brightness pattern on 
the surface of a spinning neutron star. This was con- 
firmed when burst oscillations were detected from the 
AMP SAX J1808.4-3658 close to the known neutron star 
spin frequency [1031 ]. Later, these oscillations were ob- 
served ^Jroinother AMPs near the known stellar spin 
rates jlOOl . Il04| . These established a new technique to 
measure the neutron star spin frequency. In fact, regular 
pulsations and burst oscillations are the only methods to 
measure the spin rates of neutron stars in LMXBs, and 
together they have provided spin frequencies of about 
25 sources. Apart from measuring the stellar spin, the 
modelling of shapes and amplitudes of the phase-folded 
burst oscillation intensity profiles can be useful to con- 
strain the other neutron star parameters (includi ng th e 
radius-to-mass ratio; e.g., Bhattacharyya et al. [l05l |). 
Fig. 11 shows the burst oscillation power contours and 
the power spectrum for a neutron star LMXB. The power 
contours indicate that the oscillation frequency changes 
with time. Indeed such small evolution of frequencies 
are observed for many bursts, which suggests that the 
aymmetric brightness pattern slowly moves on the spin- 
ning neutron star surfaces. It is not yet understood what 
causes such motion, although several models exist in the 
literature (27| . The oscillations during burst rise are nor- 
mally believed to be caused by the expansion of the burn- 
ing region on the stellar surface immediately after the ig- 
nition. Therefore, the study of the evolution of the burst 
oscillation properties during burst rise can be useful to 
probe the thermonuclear flame spreading. The nature 
of the aymmetric brightness pattern that causes the os- 
cillations during burst decay is not yet well understood. 
This is because the entire stellar surface is believed to be 
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ignited at this time. The existing models i nvolv e non- 
radial globa l osc illations in the surface layers [l06j , shear 
oscillations \lOl\ . etc., each of which has shortcomings. 

D. Quasi-periodic Oscillations and Broadband 
Power Spectrum 

Neutron star LMXBs show a variety of observationally 
robust timing features, many of which are not well under- 
stood yet. Continuum power, in excess of the mean ran- 
dom power, is often observed at frequencies less than 100 
Hz. The nature of this power depends on the source state 
determined from the position on the color-color plot. For 
example, a power-law- like continuum (oc v~ a ] v is the 
frequency in the power spectrum; and a is a constant) 
below 1 Hz is typically dominant for the atoll sources 
in the end part of the banana state [95[. Apart from 
the broad continua, neutron star LMXBs often exhibit 
peaked features with a wide range of frequencies (mil- 
lihertz to kilohertz). Since they are somewhat broad, 
and hence not quite periodic, these features are known as 
quasi-periodic oscillations (QPOs). Similar to the contin- 
uum features, QPOs are correlated with the color-color 
diagram. For example, a type of 15 — 60 Hz QPO appears 
in the horizontal branch of the color-color diagram of Z 
sources, and hence known as the "horizontal branch oscil- 
lation (HBO)" . The QPOs, when modelled correctly, are 
expected to significantly enhance our understanding of 
neutron star LMXBs. The observed correlations among 
different types of QPOs can be helpful for their model- 
ing. However, in this short review we will not discuss 
the QPOs and their correlations in detail, and we refer 
to an extensive review by van der Klis [95j that gives a 
large amount of information about the fast X-ray tim- 
ing properties of neutron star LMXBs. The only QPOs, 
that we would like to briefly discuss here, are called kilo- 
hertz (kHz) QPOs. These high frequency QPOs are ob- 
served from both Z and atoll sources. The frequencies 
of kHz QPOs vary between w 200 - 1200 Hz, and of- 
ten they appear in a pair (see Fig. 12). For a given 
source, the lower and upper kHz QPOs are seen to move 
up and down in frequency together without changing the 
frequency separation much. This frequency separation is 
normally found to be within 20% of the neutron star spin 
frequency, or half of that, depending on source. These 
kHz QPOs are believed to be originated from a region 
close to the neutron star, and hence, when their origin is 
fully understood, will be very useful to measure the neu- 
tron star parameters, and to test the t heory of general 
relativity (e.g., Miller et al. 1998 (lQ8j ; Stella & Vietri 
1998 [Tog; Barret et al. 2005 [ltcj). 

VII. FUTURE PROSPECTS 

Neutron star LMXBs provide a unique opportunity to 
study several aspects of extreme physics that we have 



briefly discussed in this short review. These include (1) 
constraining the equation of state models of the supranu- 
clear density matter of neutron star cores from the mea- 
surements of the stellar parameters; (2) testing the the- 
ory of general relativity; (3) understanding the physics of 
neutron star atmospheres; and (4) probing the accretion 
flow in the strong gravity region near the neutron stars. 
In addition, the X-ray study of neutron star LMXBs al- 
lows us (1) to understand the evolution of these binary 
systems; (2) to probe the long term accretion process and 
the structures of the accreted matter; and (3) to find out 
how and why these systems are different from black hole 
LMXBs. Since the observed X-ray properties mentioned 
in this review can be used as tools to achieve these scien- 
tific goals, we need to understand these properties better, 
and explore (1) the correlations among them, and (2) 
their correlations with the observations in other wave- 
lengths. These can be achieved with more sensitive X- 
ray instruments (typically with larger photon collecting 
areas), and simultaneous observations in multiple wave- 
lengths. It will also be required to monitor the entire 
sky frequently with high sensitivity in order to catch the 
transient neutron star LMXBs during their outbursts. 
Here we give a short description of some of the proposed 
future space missions with these capabilities. 



The proposed Indian multiwavelength astronomy 
space mission Astrosat will observe simultaneously in a 
wide energy range (optical to hard X-ray of 100 keV). 
Apart from this unprecedented capability, its LAXPC in- 
strument should be able to detect and measure the high 
frequency timing features, such as kHz QPOs and regu- 
lar pulsations, for the first time in hard X-rays (say, up 
to w 50 keV). Note that in its own time, LAXPC will 
probably be the only instrument capable of detecting the 
high frequency features in the X-ray band. 



The proposed next generation joint American, Euro- 
pean and Japanese satellite International X-ray Observa- 
tory (IXO) will have several highly sensitive instruments 
that will take the X-ray study of neutron star LMXBs 
to a higher level. In Fig. 13, we demonstrate how sig- 
nificant an improvement the high timing resolution spec- 
trometer (HTRS) instrument of IXO will make over the 
RXTE PCA for the high frequency timing observations. 
Moreover, Ray et al. [l 1 ll j mention that the proposed Ad- 
vanced X-ray Timing Array (AXTAR) satellite will have 
better capabilities than IXO for these timing observa- 
tions. Finally, the Japanese mission Monitor of All-sky 
X-ray Image (MAXI) will monitor the X-ray sky in the 
0.5 — 30 keV range with high sensitivity. Therefore, the 
study of the neutron star LMXBs is likely to remain an 
active and promising field in order to explore some as- 
pects of extreme physics. 
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FIG. 7: The total background-subtracted XMM-Newton RGS energy spectra (intensity vs. wavelength) of the neutron star 
LMXB EXO 0748-676 during 28 thermonuclear X-ray bursts. The data are plotted as the black histograms with la error bars. 
The blue line is the empirical continuum, with additional O VII intercombination line emission, modulated by absorption in 
photoionized circumstellar material. The upper panel is for the early burst phases (when the burst intensity is more), while the 
lower panel is for the late burst phases (when the burst intensity is less). The two red labels indicate the iron absorption lines 
at a gravitational redshift z = 0.35 plausibly originated from the neutron star atmosphere (figure courtesy: Mariano Mendez; 
Cottam et al. [H ). 
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FIG. 8: Typical color-color diagram (see text for the defini- 
tion) of the members of the "Z" subdivision of neutron star 
LMXBs. Such a diagram shows the evolution of the source 
spec tral state (figure courtesy: Peter Jonker; Jonker et al. 

El)- 




FIG. 9: Typical color-color diagram (see text for the defini- 
tion) of the members of the "atoll" subdivision of neutron star 
LM XBs (figure courtesy: D iego Altamirano; van Straaten et 
al. [l!4| : Altamirano et al. [1151 ]'). 
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FIG. 10: Power spectrum (i.e., Fourier transform of the observed X-ray intensity profile) of a neutron star LMXB. The peak 
at 442.36 Hz implies that (1) there is a periodic intensity variation with that frequency; (2) this source is an accretion-powered 
millisecond pulsar (see text); and (3) the spin frequency of the neutron star is 442.36 Hz (figure courtesy: Diego Altamirano; 
Altamirano et al. [l!6l ]). 



20 



K 
o 

CJ 

o 



533 
532 

531 

530 

529 

528 

527 
526 




2 



6 8 
Time (s) 



10 12 



o 

4 8 



- 2 



o 




10 100 1000 

Frequency (Hz) 



FIG. 11: Upper panel: X-ray intensity profile of a thermonu- 
clear X-ray burst from the neutron star LMXB 1A 1744-361 
[l!2j |. Power contours around 530 Hz (calculated by Fourier 
transform) implies that there is a significant intensity varia- 
tion (burst oscillation) with « 530 Hz frequency. Note that 
the shape of the power contours indicates a slight change in 
frequency with time, which is common for burst oscillation. 
Lower panel: The corresponding power spectrum. The peak 
at ~ 530 Hz implies the burst oscillation with that frequency. 
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FIG. 12: Power spectrum (i.e., Fourier transform of the ob- 
served X-ray intensity profile) of a neutron star LMXB show- 
ing a pair of kHz QPOs (see text). The vertical dotted lines 
show the centroid frequencies of these QPOs. 
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FIG. 13: A comparison of a power spectrum observed with 
two (one current and one future) instruments. The lower 
panel shows that the lower kHz QPO is barely detected and 
the upper kHz QPO is not detected from the power spectrum 
created from an intensity profile of 100 seconds duration ob- 
served with the PCA instrument of the RXTE satellite. The 
upper panel shows that, if the same observation were made 
with the High Time Resolution Spectrometer (previously pro- 
posed for the X-Ray Evolving Universe Spectrometer (XEUS), 
but now proposed for the International X-ray Observatory 
(IXO)), both the kHz QPOs would be very significantly de- 
tected (figure courtesy: Didier Barret; Barret et al. [ll7j ]). 



